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Advances in the field of numerical relativity now make it possible to calculate the final, most 
powerful merger phase of binary black-hole coalescence for generic binaries. The state of the art has 
advanced well beyond the equal-mass case into the unequal-mass and spinning regions of parameter 
space. We present a study of the nonspinning portion of parameter space, primarily using an 
analytic waveform model tuned to available numerical data, with an emphasis on observational 
implications. We investigate the impact of varied mass ratio on merger signal-to-noise ratios (SNRs) 
for several detectors, and compare our results with expectations from the test-mass limit. We note 
a striking similarity of the waveform phasing of the merger waveform across the available mass 
ratios. Motivated by this, we calculate the match between our 1:1 (equal mass) and 4:1 mass-ratio 
waveforms during the merger as a function of location on the source sky, using a new formalism for 
the match that accounts for higher harmonics. This is an indicator of the amount of degeneracy in 
mass ratio for mergers of moderate-mass-ratio systems. 

PACS numbers: 04.30.Db, 04.80.Nn 95.30.Sf, 97.60.Lf 



I. INTRODUCTION 

The merger of a black-hole binary will be one of the 
strongest sources of gravitational waves, with a greater 
luminosity than the combined electromagnetic luminosity 
from all the stars in the visible universe. Ground-based 
detectors like LIGO, Virgo, and GEO, currently entering 
their second generation of development, are sensitive to 
the mergers of stellar black holes, while the space-based 
LISA will observe mergers of massive and supermassive 
black holes. It has long been expected that the final 
mergers of black-hole binaries would be significant for 
interpreting gravitational wave measurements. While a 
physically motivated model, and corresponding template 
bank, would not be necessary for detection such a 
model would be the only avenue toward extracting all of 
the available information about the system that is con- 
tained in the merger signal. And since the merger is 
likely to constitute the majority of the detectable sig- 
nal for the next generation of ground-based detectors , 
such a physically motivated model would be necessary 
for gaining an understanding of the physical sources gen- 
erating the detected signals. 

In the absence of merger models, early investigations 
had to use information from perturbative approximations 
to guess at the impact of mergers. In the New- 
tonian approximation for the gradual adiabatic inspiral 
of the holes, combined with the understanding of the 
post-merger ringdown as the quasi-normal modes of a 
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Kerr black hole, was used to guess at the contribution of 
mergers to the signal detectability. This guess was essen- 
tially validated by the observed behavior of numerically 
simulated merger signals [2[. However, while the power 
spectrum could be approximated, the physics behind the 
power spectrum, the amplitude and phase evolution that 
would lead to that spectrum, the accuracy with which 
the merger phase could be simulated or modeled, and 
the amount of information about the source that could 
be extracted from detected signals were completely open 
questions. Over the course of the last few years, the field 
of numerical relativity has provided a means of studying 
the detailed structure of these merger signals for the first 
time. Initially focusing only on the equal-mass, nonspin- 
ning case, several groups have since explored both the 
nonspinning axis of parameter space as well as the vast 
expanse of spinning parameter space 04^- The current 
availability of merger waveforms now makes it possible 
to address the questions previously mentioned, through 
the measurement of these signals. 

Because the merger is the dominant contributor to the 
overall signal power, particularly for ground-based de- 
tectors where it provides the majority of the detectable 
signal, answering these questions is a critical exercise ac- 
tively being addressed by many groups. Significant atten- 
tion has been given to the problem of modeling the sig- 
nals with sufficient accuracy for detection with ground- 
based observations. Recent work has begun addressing 
not only detection, but also estimating the source param- 
eters using ground-based [l^l and space-based detectors 
(Tll - ITi) . Much work has also gone into interpreting the 
available merger waveforms, in an attempt to better un- 
derstand them, both with regard to a physical intcrpreta- 
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tion of the source [15l | and with regard to understanding 
what drives the recoil of systems due to asymmetric ra- 
diation [iMl. 

In this paper, we revisit the nonspinning subset of pa- 
rameter space, with the goal of studying the observa- 
tional implications of nonspinning merger waveforms. In 
Sec. ini we briefly describe the procedure for generating 
complete nonspinning waveforms. In Sec. IIIH we study 
the contribution to the achievable SNR from the inclu- 
sion of mergers, and its variation with mass ratio. In 
Sec. IIVI we study more detailed comparisons of the non- 
spinning waveforms, including a novel implementation of 
the "match" statistic [ll|. In Sec. |Vl we discuss the 
general observational implications of nonspinning merger 
waveforms. In the Appendix, we derive the formalism for 
the novel match implementation employed in Sec. IIVI 



II. METHODOLOGY 

The observable quantity being measured by gravita- 
tional wave interferometers, be they ground- or space- 
based, is the strain on the spacetime, h = 6L/L, or its 
derivatives. We therefore require a model of the wave- 
form, h(t), that we expect to measure from black- hole 
binaries. The models we employ are predictions for the 
emitted strain, or the strain in the source frame. The 
detected strain depends on the distance to the source, 
the position on the detector's sky, and the detector's re- 
sponse. 

In this work, we focus on Advanced LIGO for ground- 
based observation, and LISA for space-based observation. 
For Advanced LIGO, we assume a constant response as 
a function of frequency, which should be adequate for all 
but the lowest-mass cases. For the detector noise, we use 
the wide-band tuning [201 typically associated with burst 
sources as was done in [2|, due to the superior sensitiv- 
ity at higher frequencies where the merger will occur for 
lower masses. For LISA, we employ the effective noise 
floor from [2l|, [22| , which includes both the contributions 
from noise sources as well as the average response of the 
detector to signals and instrumental noise, which is non- 
trivial for the higher frequencies in LISA's band. For 
frequencies in the range 3 x 10~^IIz < / < 1 x 10~"'llz, we 
employ a more conservative estimate of the acceleration 
noise as was done in 0, instead assuming a steeper am- 
plitude spectral density that falls off as Below 
3 X 10~^Hz, we assume the detector has no sensitivity. 
We apply an overall factor of 3/20 to the LISA power 
spectral density as discussed in [2J]. 

The remaining element is the model of the emitted 
waveform, which will depend on the intrinsic parameters 
(i. e. the mass, mass ratio, spin vectors, and eccentricity), 
and which will vary over the sky of the source. The 
emitted waveform can be conveniently represented by a 
harmonic mode decomposition. If h is the complex strain. 



then the mode decomposition is given by 

oo e 

h^Y^Y. h,ra{t,R)-^Ytm{e,4>), (1) 

where ~^Yfm are the spherical harmonics of spin- weight 
(-2) [2^. Being complex, h contains both wave polar- 
izations, defined by the relationship h = + ih^- For 
an equal-mass system, h{t, R, 0, (p) is dominated by the 
quadrupole, the combination of^ = 2,TO = +2 modes: 

/iquad = h22 -^Y22{e, 0) + /l2-2 "'i^2-2(^, ■ (2) 

Additionally, symmetry considerations for equal-mass 
nonspinning systems demand that /122 = ^2-2! therefore 
wc will often use /122 as a proxy for the full quadrupole 
waveform. 

For the waveform comparisons presented in this work, 
we use the model first presented in [l^, which has been 
validated by comparison with available data from nu- 
merical simulations for all harmonic components through 
£ = A. The model, referred to as the IRS-EOB model, 
uses the efFcctive-onc-body (EOB) Hamiltonian formal- 
ism for the inspiral [26j . For the merger-ringdown, we 
employ a novel paradigm which we call the implicit rotat- 
ing source (IRS), wherein we apply a fit to a physically- 
motivated functional form for the phasing (Eq. 9 in [l^ ) , 
and, for the amplitude, a model for the fiux constrained 
to be consistent with the inspiral flux through 3.5 post- 
Newtonian (FN) order and to vanish as it approaches the 
ringdown frequency (referred to as "Model 2" and given 
by Eq. 19 in [15|). In Figure [H wc compare the intrinsic 
error in phase and amplitude for the model and for our 
numerical waveforms, using 4:1 as a representative case. 
For brevity, we will refer to the unequal-mass runs as ra- 
tios, i. e. the q = M1/M2 = 1/4 run will be the 4:1 run, 
and the ratio notation will only be used in this context. 
We shift all waveforms to peak in amplitude at t = 0, and 
to agree in phase a.t t — — 500A/. Wc use the difference 
between our two highest resolutions as an indicator of 
our numerical error, and we assume errors in the model 
due to different model parameters arc independent, and 

calculate a phase error 5(j) ~ (^'§X'^-^iJ ' where Xi 

is simply the A parameter as used in [2^1 for the EOB 
inspiral, and = {«;, 6, tp, Clg, flf} for the IRS merger, 
using the notation in [l5|. For 5Xi, we use the values in 
Table II of dl. 

Figure [2] shows the quadrupole radiation for four mass 
ratios - 1:1, 4:1, 6:1, and 20:1 - generated using the IRS- 
EOB model. Wc note that in using the 20:1 mass ra- 
tio, we have extrapolated to mass ratios that cannot, as 
yet, be validated by simulation. The amplitudes for all 
the runs in Figure [D have been rescaled to better agree 
with the equal-mass amplitude, using the leading-order 
Newtonian scaling. This also emphasizes the phasing 
agreement that begins in the late inspiral and contin- 
ues through the merger waveform, which was discussed 
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FIG. 1: Comparison of amplitude and phase errors between 
the numerical data and the IRS-EOB model, for the case of a 
4:1 waveform. The inherent phase inaccuracy of the model is 
significantly smaller than the numerical phase error, but the 
amplitude errors are comparable for the model and numerical 
data. 



in [15| and will be the topic of further discussion in a 
later section. 

Since many investigations relating to both LIGO and 
LISA have focused on detectability, rather than charac- 
terizing the signal, a model of a quadrupole-only signal 
has been adequate within their margin of error. However, 
higher harmonics can be more significant for calculations 
such as determining template fidelity with the match, or 
any attempt to extract source parameters from the sig- 
nal, as such investigations depend sensitively on the fine 
detail of the phase evolution. We will investigate the im- 
pact of higher harmonics in this context, and we will also 
include higher harmonics in our calculations of SNR, al- 
though the SNR contribution is essentially negligible for 
all cases investigated here with the possible exception of 
20:1. 



III. SNR AND POWER SCALING 

SNR is the most useful statistic for assessing the de- 
tectability of a given signal with a particular detector. 
The SNR, which we denote as p, is given by 

- {h\h), (3) 

where "(-I-)" denotes a noise-weighted inner product, 
given by 



{h,\h2}=2 J df 
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and Sn{f) is the power spectral density of the detector 
noise discussed in the previous section. The sky-averaged 
SNR is given by 



{p')= jd{\nf) 
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Here, /ichar(/) = 2/|/iopt(/)| is the characteristic sig- 
nal strain, and = A/5/irms(/) = \/5/5„(/) is the 
root-mean-square of the detector noise fluctuations mul- 
tiplied by \/b for sky-averaging. /iopt(/) is the Fourier 
transform of the optimally-oriented signal strain [ij . For 
quadrupole-only cases, the orientation-averaged signal 
strain is trivially calculated from the optimally-oriented 
strain by dividing by y/Z. 

Before the advent of merger waveforms from numerical 
relativity, expectations about the power scaling of the 
merger waveforms, and thus the achievable SNRs, were 
formed by using the test-mass limit as a surrogate, while 
the scaling of the inspiral power can be approximated 
by FN expansions in the weak-field limit. Specifically, 
we know that the SNR from the inspiral scales as ^ to 
leading order, where 77 = MiM2/M'^ is the symmetric 
mass ratio of the binary. It was further assumed in 
based on the prediction for total radiated energy in the 
test- mass limit that the merger SNR scales as ry. We 
note, however, that in [31, the peak (2,2)-mode energy 
flux was best fit by the function 



E22 = 4.40 X lO^^Ty^ 



5.43 X lO^^T]^ . 



(6) 



with £^22 oc |/i22p oc p^. The signiflcant improvement 
in performance of the 77^ + rj'^ fit, compared to a strictly 
1]^ fit, may indicate that differences between the physics 
of the merger for comparable masses and the test-mass 
plunge are being measured. The absence of a well-defined 
innermost stable circular orbit (ISCO) in the equal-mass 
case [2^ [13, compared to the obvious ISCO for suffi- 
ciently small mass ratios, further supports this picture. 
Certainly as the masses differentiate more, the test-mass 
analogy bears out more. For the equal-mass case, Eq. ([6]) 
indicates that the 77^ and rj'^ terms contribute roughly 
equally, but the 77^ term obviously becomes less impor- 
tant for ever-smaller mass ratios. 

The different scalings with rj are illustrated in Figure |31 
where we plot the Fourier transform of hybrid waveforms, 
constructed analogously to [2| by tying a FN inspiral to 
our numerical data at a point where they reach equal ac- 
curacy. We do this in part to emphasize that the change 
in scaling of the merger signal is not an artifact of the 
IRS-EOB model, but is apparent in our raw numerical 
waveforms. The left panel demonstrates the ^/rj scaling 
of the inspiral (Mw < 0.08) for the 1:1, 2:1, 4:1, and 
6:1 quadrupolar waveforms, and the right panel shows 
that the merger scaling is well approximated by a lin- 
ear dependence on rj for the merger (Mw > 0.08). The 
deviation of the peak E22 from a simple quadratic-in-77 
scaling appears to be due primarily to differences in the 




FIG. 2: In the top panel (a.), we show quadrupole waveforms generated by using the model presented in [isl ] for mass ratios of 
1:1 (equal mass), 4:1, 6:1, and 20:1. When the waveforms are aligned in time based on their peak amplitudes, and aligned in 
phase to agree at said time, there is significant overlap of the waveforms for the 1:1, 4:1, and 6:1 cases over the final ~ 5 cycles 
leading up to merger, which is shown more clearly in the bottom panel (b.). 



frequency of the peak, since the signals closely follow a 
linear-in-77 scaling when evaluated at the same frequency 
prior to the peak. 

For equal masses and moderate mass ratios, 
h{t,R,9,(f>) « hquad{t, R), so that averaging the SNR 
over the binary's orientation is trivial. This may be sufB- 
cient for ground-based detectors if they are primarily de- 
tecting stellar-mass black holes {M < IOOMq) due to the 
limited available mass range. If intermediate-mass black 
holes (IOOMq < M < IO^Mq) exist, then smaller mass ra- 
tios {q <C 1) may occur, and higher-order modes will con- 
tain progressively more power relative to the quadrupole 
as (? — > 0. Rather than ignoring higher-order modes or 
introducing complexity by averaging over them, we in- 
stead focus on the optimal orientation of the binary, and 
only average over the sky of the detector. 

To demonstrate the decrease in SNR with a signifi- 
cant deviation from equal mass, we show the 1:1 and 
6:1 cases for Advanced LIGO and the 1:1 and 20:1 cases 
for LISA in Figures 0] and [51 respectively. The panels 
show contour lines for both mass ratios, with one set of 
lines corresponding to the SNR accumulated before the 
corresponding Schwarzschild ISCO, and the other cor- 
responding to the full signal. As described earlier, the 
SNR decreases as rj deviates from 0.25, with the inspiral 
SNR scaling as ^/fj and the merger SNR scaling as 77. We 
therefore expect the SNRs to scale roughly as y/rj for the 
lowest masses where the inspiral matters most, and as 77 
(potentially with terms of higher power in 77 as well) for 
higher masses where the merger contributes the majority 



of SNR. 

In both mass-ratio cases for each detector, the late 
inspiral-merger phase constitutes the majority of the 
SNR in high-SNR events. The merger contributes sig- 
nificantly for masses M > SOM© for Advanced LIGO, or 
(1 -I- z)M > IO^Mq for LISA. As previously reported for 
equal-mass mergers [l|, 0] , the merger contribution to the 
signal tends to dominate strongly for these larger-mass 
systems. For the unequal-mass cases the merger plays 
a dominant role over a similar range of masses, though 
the level to which the merger dominates the overall SNR 
is significantly diminished for very unequal masses (right 
panels of Figures H] and [5|) compared to the equal-mass 
case (left panels; see also the figures presented in [H, 0]). 
In some ways the equal-mass case is exceptional, rather 
than representative. For observations of IMBH merg- 
ers {M > lOOAf©) with Advanced LIGO, however, the 
merger always dominates, as the relatively sharp wall in 
low-frequency sensitivity effectively wipes out the inspi- 
ral contribution. 



IV. WAVEFORM COMPARISON 

Numerical relativity now provides a reasonably clear 
picture of the late stages of merger, specifically in the 
form of the waveforms. While considerable progress has 
been made in understanding how to detect inspiral sig- 
nals and characterize how they depend on system param- 
eters, including mass ratio and spin, there is little simi- 
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FIG. 3: Scaling of the Fourier amplitude of hybrid (PN inspiral/numerical merger) waveforms for different mass ratios by yfr\ 
in the left panel (a.), and by t] in the right panel (b.), which appears to be an excellent approximation. 




M{Mq) M{Mq) 

FIG. 4: SNR contours for Advanced LIGO with g = 1 (a.) and q = 1/& (b.). Note, when comparing the two panels, that the 
masses are total masses, which determine the overall waveform amplitude. The solid lines correspond to the SNR calculated 
from the full waveform, including the merger, while the dotted lines correspond to the SNR contribution from the portion of 
the signal with frequency lower than the Schwarzschild ISCO frequency. 



lar work addressing observations over the signal-space of 
mergers. Most observational work so far has considered 
these effectively as unmodeled sources. In [l^, we exam- 
ined the relationships between the merger waveforms and 
the physical motion of the source, emphasizing simple 
common features in order to form a general characteri- 
zation of nonspinning mergers. These features have ob- 
servational consequences as well. In particular, we noted 
general similarity in the late-time portions of dynamics 
and waveforms, which we now revisit. 

In Figure |6] we plot model waveforms for the 1:1 case, 
the 4:1 case with its amplitude rescaled with the leading- 
order T] dependence, and the difference 5h = hi — /12 
between these two waveforms. By doing so, we see that 
the apparent phase agreement shown in [15| is, not sur- 
prisingly, partially an artifact of aligning all the phases at 
the peak strain amplitude, and thereby enforcing a node 



in 5h at that time. However, the merger is unique in the 
suppression of the final beat prior to ringdown in 5h, and 
enhancement of that beat if a 7r/2 phase shift is applied. 
This can be seen in the time series of Figure |6l but is 
most evident in the Fourier representation of Figure [T] 
This extended frequency agreement also provides a sim- 
ple explanation for the ^Jrj difference in scaling between 
the inspiral and merger. The amplitude scales linearly 
with to a good approximation for both the inspiral 
and merger, but the time interval spent within a given 
frequency bin scales as rj~^/^ to leading order for the 
inspiral and is nearly constant for the merger for moder- 
ate mass ratios. This results in a relative 77^/* amplitude 
scaling of h, and therefore SNR, between the inspiral and 
merger. 

While SNR is certainly the most relevant statistic for 
detection purposes, it tells us little about the details of 
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FIG. 5: SNR contours for LISA with q = 1 (a.) and q — 1/20 (b.). The sohd hnes again correspond to the full waveform 
SNR, while the dotted lines correspond to the SNR contribution from frequencies lower than the Schwarzschild ISCO frequency. 
While the observable range of high-SNR mergers is reduced by a factor of several at 20:1 from what was seen in the 1:1 case, 
sources are still easily detectable to large redshift over a similarly broad mass range. 




FIG. 6: Differences, Sh, in 1:1 and rescaled 4:1 mass ratios. The inspiral waveform evolves out of phase on a timescale shown 
clearly by the beats in 5h (a.). The phase alignment persists for a slightly longer time during the merger (b.). 



the waveform, in particular the evolution of the phase, 
which may be critical when answering questions regard- 
ing signal characterization. The "match" jl^ is a useful 
statistic for more detailed waveform comparisons, as it is 
sensitive to small differences in waveform phase. For any 
two waveforms, hi and the match A4 is defined using 
the noise- weighted inner product 

M = ^±^^=. (7) 

^J{hi\hi){h2\h2) 

The match can be viewed as the fraction of the matched- 



filter SNR that is recovered by using /12 as a filter to 
search for hi, rather than using hi itself (the optimal 
filter) . The left panel of Figure |S] shows a typical com- 
parison for the 6:1 case, which should have the strongest 
higher harmonics among the numerical simulations stud- 
ied here. Also, we show a frequency-based comparison 
for the same case in the right panel of Figure [51 where 
the (2, +2) modes can be seen to dominate the signal 
power until well into ringdown (indicated here by a ver- 
tical dashed line). Nonetheless, it is still possible that a 
sub-dominant mode may modulate the signal to a suf- 
ficient degree to significantly diminish the recoverable 
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FIG. 7: The Fourier amplitude of the waveform difference in 
Figure [S] both with the phase shift shown there, as well as 
shifted by n/2 to illustrate the dependence of the power at 
high frequencies on the phase alignment. 



SNR, or to impact the template member that has the 
highest likelihood, if only the dominant mode is used as 
a filter. We will therefore develop an appropriate formal- 
ism for including all modes analytically in a calculation 
of the match. 

We can further calculate the SNR of the difference in 
waveforms, 5h (sec Figure IH]), which is essentially a mea- 
surement of our ability to distinguish two waveforms from 
each other. This simple statistic is related to the "mis- 
match" , 1 — , as well as the SNR, p (see also [s^] ) : 

{Sh\6h) = {hi - h2\hi - h-i) 

= {hi\hi) + {h2\h2)-2{hi\h2) 

« 2p'{l-M), (8) 

where \hi\ = ^ lji\\h\)^ and the final approximation 
comes from assuming that the SNRs oi h\ and are 
approximately equal. Wc note that the curves "1:1", 
"rescaled 4:1" and "1:1 - 4:1" in Figure [ 71 arc simply 
the integrands of |/i2| and \J {dh\Sh), respectively, 
without noise- weighting. The latter is roughly an order 
of magnitude smaller than \hi \ during the entire merger 
phase. In this case we can see from Eq. (|8]) that the match 
among moderate-mass-ratio mergers is likely to be quite 
high for ground-based interferometers. We can there- 
fore expect that, for instance, a small subset of merger 
waveforms would be capable of sufficiently covering a 
large range of nonspinning parameter space for detec- 
tion purposes, but the apparent mass-ratio degeneracy 
in the merger will have a negative impact on parameter- 
estimation efforts. 

Thc SNR has a trivial inverse proportionality with lu- 
minosity distance, so swapping SNR and luminosity dis- 
tance can give you, for instance, the distance horizon at 



which the difference between waveforms can be detected, 
by setting the SNR at some fixed threshold. We do so 
in Figure ini using the same 6h as above as an example. 
We use a fixed SNR of p = 10 as the threshold of de- 
tectability. The values in Figure [9] have an interesting 
implication in light of Eq. , in that for sources farther 
than the distance horizon, we cannot distinguish between 
a 1:1 waveform at Dl, or a waveform with mass ratio rj 
at a distance {i]/0.25)Dl. 



To further investigate the implications of the appar- 
ent degeneracy in mass ratio for moderate-mass-ratio 
mergers, we calculate the match between the full 1:1 
and 4:1 waveforms, including all available harmonics, as 
a function of the source orientation. In Figure I10[ we 
show as an example a comparison with redshifted mass 
(1 + z)M = 3 X 10*^Afo for LISA and M = lOOMo for 
Advanced LIGO, although the result will be qualitatively 
similar for any masses where thc merger is emphasized 
relative to the inspiral in the whitened waveform (see the 
Appendix and Figure [TT|) for a given detector (i. e. any 
case where the signal merges at or below the peak sen- 
sitivity of the detector) . Figure [10] shows a sky map of 
the matches for Advanced LIGO and for LISA, where 
the match at each point corresponds to the equal-mass 
waveform calculated at that point on the sourcc sky. The 
match is maximized over the orientation of the 4:1 wave- 
form in the (f> direction at a given 9, where </> is the az- 
imuthal coordinate and 8 is the polar coordinate. The 
maximum at a given 9 then corresponds to the "best" 
match, and the minimum to the "minimax" match (3l| . 
The azimuthal sky position is degenerate with the orbital 
phase, so that the maximization procedure is identical to 
finding the maximum value in the azimuthal direction for 
a particular inclination. We do not maximize over the in- 
clination, since the spin- weighted spherical harmonics are 
a more complicated function of polar angle. In the sky 
maps, it is clear that the maximization over polar angle 
would occur at thc poles, where the quadrupole modes 
are most dominant. This is consistent with our previous 
results showing the striking similarity of the quadrupole 
radiation across modest mass ratios. We observe the ex- 
pected "north/south" symmetry, since all nonspinning 
binaries evolve in a fixed plane. The azimuthal asymme- 
try is greatest in the orbital plane, where the fractional 
luminosity of the higher harmonics relative to the dom- 
inant quadrupole modes is greatest. We note that the 
sky map would be uniform for single mode matches, so 
the structure in Figure [TUj is the result of the harmonic 
content, and therefore requires the formalism contained 
in the Appendix in order to maximize quasi-analytically. 
The average match over the sky of the source for the 
cases in Figure [TUj is 0.96 for Advanced LIGO, and 0.95 
for LISA. This means that the 1:1 waveform can be con- 
sidered an effective template (in the sense of [3l|) for 
typical Advanced LIGO mass ratios for a large fraction 
of sourcc orientations. 
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FIG. 8: In the left panel (a.), we show a time series representation of a 6:1 mass-ratio 10^ Mq black-hole binary at a distance of 10 
Gpc, using the waveforms from our numerical simulation. This is a typical example of the (2, + 2) modes constituting the vast 
majority of the overall power content of the waveform. This is further demonstrated in the Fourier-series representation shown 
in the right panel (b.), where the (2, +2) modes dominate until well into the ringdown, the onset of which is approximately 
indicated by the dashed vertical line. 




FIG. 9: Luminosity distance where the difference between the 1:1 and rescaled 4:1 mass ratio waveforms is detectable with an 
SNR of 10 for Advanced LIGO (a.) and LISA (b.), also referred to as the distance horizon. This can be interpreted as being 
the maximum distance at which we can distinguish these two sources with each interferometer. 



V. CONCLUSIONS 

We have applied a model for nonspinning late inspiral- 
merger-ringdown waveforms to answer questions regard- 
ing the implications of including the merger phase in data 
analysis efforts. We have verified that, while the merger 
contributes a smaller fraction of the total SNR as we 
deviate from the equal-mass case, it still dominates for 
moderate mass ratios, providing nearly the entirety of the 
detectable signal for ground-based observations of IMBH 
systems. In addition, we have studied the commonality 
previously observed in the phase evolution of the merger 
waveform for moderate mass ratios. While this com- 
monality bodes well for detection, since the equal-mass 
merger waveform alone would do well as a search filter 



for all moderate mass ratios, this has negative implica- 
tions for signal characterization. Indeed, by calculating 
the "match" as a function of location on the sky of the 
source, we have demonstrated that the equal-mass wave- 
form can be considered an effective template for detect- 
ing other moderate mass ratio signals for a wide range of 
source orientations for both Advanced LIGO and LISA. 
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FIG. 10: Sky map of matches between the 1:1 (equal mass) and 4:1 waveforms for Advanced LIGO (a.) and LISA (b.). 
Throughout the text, if is the azimuthal coordinate and 6 is the polar coordinate. The specific case shown corresponds to a 
mass of M — IOOMq for Advanced LIGO, and (1 + z)M = 3 x 10^ Mq for LISA (the luminosity distance is irrelevant for this 
calculation), but the result will be qualitatively similar for any mass M > IOOMq for Advanced LIGO and M >2x 10® Mq 
for LISA, due to the constancy of the detector response and the similar spectral content of the noise over the band of the 
signal for those cases, i. e. cases where the "whitening" procedure (see the Appendix) emphasizes the merger. The sky location 
corresponds to the position on the sky of the 1:1 waveform, with the 4:1 waveform being rotated in the 4> direction to maximize 
the match. Therefore, for a fixed 9, the maximum in will correspond to the "best" match, and the minimum to the "minimax" 
match. 



atcd Universities. 



Appendix: Generalized Phase Maximization 

We are interested in generalizing the procedure pre- 
sented in jsH for maximizing the match ([7]) vifith respect 
to the initial orbital phase constants between a target or 
exact (label X) waveform and a template or approximate 



(label A) waveform. Specifically, whereas the previous 
method is restricted in its validity to radiation that is 
quadrupole-only, we wish to derive the general method 
for maximizing the match for arbitrary harmonic content. 
Wherever possible, we preserve the original notation from 




For the exact and approximate waveforms, we can rep- 
resent the measured strain waveform as 
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= SR 



3? 



= -P'l'^y^mll^fai^l [cos(m0) cos(TOiy9'^''^) + sin(m0) sin(m(p'*''^)] 



(9) 



where A;^^ ee cos(mv3^''^), = sin(m(^^''^), /ii£„i = 

F\hira\cos{m(l)), h2im = F| | sin(m(/)) , Fe'" = F+ + 
ii^x is the complex beam pattern function, and 6 and (p 
describe the angular position on the source's sky (with 
K absorbed into the definition of <j)). Since (p can be ab- 
sorbed into V?, the following procedure maximizes over 
the relative azimuthal orientation as well as the orbital 
phase. For this analysis, we assume a common source 
polar angle 9 for the exact and approximate waveforms, 
although the procedure could be further generalized to al- 
low maximization/minimization over all relevant angles. 

As in [3l[, we wish to find the "best" and "minimax" 
match. We therefore wish to form an appropriate basis 
in which we can decompose the exact and approximate 
waveforms separately, and subsequently find the projec- 
tion of the resulting approximate "vector" on the exact 
"vector". Conceptually, in [3l| the procedure amounted 
to finding the ellipse resulting from projecting the circle 
that the approximate waveform makes in its 2-plane onto 
the 2-plane formed from the decomposition of the exact 
waveform, where the 2-planes are the spaces spanned by 
the orthonormal bases constructed using the exact and 
the approximate waveforms. In our analysis, we extend 
this concept to finding the minimum and maximum ra- 
dius resulting from finding the sum of projections of ap- 
proximate circles for the available modes on the exact 
planes corresponding to those same modes. We could 
alternatively include in the sum the cross contributions 
from particular approximate circles for a given mode on 
the planes of all available exact modes. However, while 
such a result would be more directly related to the SNR 
achievable by using the approximate waveform as a tem- 
plate, the inclusion of cross-mode contributions would 
be unphysical and less useful as a gauge for potential pa- 
rameter estimation. We therefore include only like-mode 
contributions, although the following derivation can be 
trivially altered to include all cross-mode contributions, 
and the final result will be the same in all but the most 
exotic cases. 

To form the desired bases, we first construct a pair of 
"whitened" vectors [sj , as shown in Figure [TI] in both 



the approximate and exact planes, to account for the 
presence of noise, the detector response to noise, and the 
detector response to the raw signal h/^'/^ (where n is 1 or 
2), 

htl= J df^e-^'-f\ (10) 

— oo 

where "h" denotes the Fourier transform of h. With these 
whitened vectors, the noise-weighted inner product (|3]) 
can be easily calculated in the time domain: 

/•oo 

(/ll|/l2) = / dth'*{t)h'2{t). (11) 

Jo 

Instead of attempting to construct a single orthonor- 
mal basis, we generate an orthogonal (not normal) basis 
for each mode, with the normalization chosen so that 
the sum over modes is normalized, i.e. 

lA.X' 

A,x ^ ri-um 

,A,X" _ ,A,X'_nA,X', A,X\^ A,X 

4t = I ^ ■ (12) 

This expression yields an appropriate normalization over 
all modes, since e^'-^ = X]^m normalized by con- 

struction, with the individual £m modes being appropri- 
ately weighted by their relative barycentric power and 
the response of the detector. Eq. p2|) is therefore a set 
of orthogonal bases which are all constrained by the to- 
tal signal power, and by their common dependence on the 
orbital phases of the exact and approximate waveforms. 
We therefore retain the original two degrees of freedom 
as in [sij . 

If we focus only on like-mode contributions to 
the match, we can construct the projection operator. 
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FIG. 11; Examples of "whitened" waveforms [s^l that we use to form a basis for calculating the match. Examples for Advanced 
LIGO are shown in the left panel (a.), and examples for LISA are shown in the right panel (b.). The ordinate values are 
arbitrarily scaled. 



Px{e-^), and the resulting projection pa, of a vector where 
onto the X-plane, 



_X 



2 e„ 



where is shorthand for , and 

nim 

defined as 



(13) 



A.X • 



n=l 1=2 



A,X I \ A,X , ■ I \ A,X /-.^x 

e^im = cos(ma) e^l^ + sin(ma) e^l^ , (14) 



where a is an arbitrary initial orbital angle for the ap- 
proximate waveform. Substituting for e^^^ from ()14| 
into ([T^ yields 



9a = ^ [cos(ma) pifon + sin(ma) p2foi 



(15) 



where 



Vnim = -Px(e„f,„) 



I X I X \ Urn ■ 
k=l \^klm\^kim/ 



(16) 



If we again focus only on like-mode contributions for sim- 
plicity, then \pa\ can be expressed in a form which rep- 
resents, geometrically, a sum of ellipses, given by 

bap = ^ [Aim cos'^ {ma) + Bini sin^(TOQ;) 



em 



+2 Cim cos{ma) sin(TOQ;)] , (17) 



A 



im 



\P 



.A 



Be 



'2im I 
A 



bwm) 



.X \2 



(e. 



X 



-'2tmV--\lml 



(e 



^x 



-■2lr, 



.X \ 
^llm/ 



(e: 

mj 

i^Uml^Umf 



.X \ 



(e 



'Um/ 



+ 



(e 



■-Um 
\^2em 



--Un 



.X 
-■2111 



^X 



.X \ 



(18) 



-'2£mP2£m/ 



While (jl7p is trivial to maximize or minimize analyti- 
cally for the ease of a single mode as in [3l| , the case of 
multiple modes generally requires a numerical solution. 
However, if we assume a single mode (or mode pair) is 
significantly larger than any other mode, then we can 
specify an approximate solution for the value of a that 
yields the "best" match. Generally, the condition for ex- 
tremizing (|17p is given by 

[m {Aim - Bim) sin(2mQ;) 

-2 m dm cos(2ma)] = . (19) 

We can then apply the aforementioned assumption that 
a single mode pair dominates. In geometric terms, this 
means that we assume that the semi-major axis for the 
dominant mode(s) in Eq. ([T7)) is larger than the quadra- 
ture sum of the semi-major axes of all other modes. In 
this case, the largest value for Eq. ([T7|) will occur very 
near the a that maximizes the dominant mode(s), with 
the other modes providing at most a small perturbation. 



12 



This condition can be expressed as 



E 



2 



im^CM 




^ CM - 



+ CL> (20) 



where CM corresponds to the dominant mode(s), and we 
only include the larger roots of Eq. ([TO)) corresponding 
to the "best" match. For all cases in this paper, C = 
\M\ = 2, with \paCM\ = \PaCi-M)\ by symmetry, so 
that the condition for ([20]) in this case will be £m ^ 
li\M.\. Finally, we can calculate the condition on a for 
maximizing the match under these assumptions: 



ttbost 



2M 



cos 



{A 



CM 



Bcm) 



^{AcM-BcMf+AC 



1 

CM , 



.(21) 



We reiterate that Eq. ((2T|) is not valid if subdominant 
modes contain comparable power to the dominant mode 
or mode pair, and a similar method cannot be used to 
find the "minimax" a. In these cases, Eq. (|T^ can only 
be solved numerically. Even if a numerical solution is 
required, this method is still more efficient than a brute- 
force maximization over (p-^ and , as it makes it a one- 
dimensional search over a. Since ^£,„, i?^m, and Cf,„ are 
all less than unity, the error in the match will be of the 
same order as the sampling interval in a over the range 
(0, 27r], assuming modes with very large m are negligible. 
For instance, in this work we only include modes. 



with m = + 4 the largest relevant m mode. Therefore, in 
order to calculate the match to three significant digits, 
we take 10^ samples of a and record the global max- 
imum and minimum, corresponding to the "best" and 
"minimax" matches, respectively. One could implement 
a more clever algorithm, such as Brent's method, if the 
required accuracy for the match makes the sampling pro- 
cedure too computationally expensive. We have verified 
that the minimum and maximum from Eq. (|19p agrees 
with the minimum and maximum found using a Nelder- 
Mead simplex over the two-dimensional ip-^-kp^ space. 

We note that, even if we were to include the cross- mode 
contributions in our derivation, the result would remain 
a sum of a set, albeit a much larger set, of ellipses. It 
should be noted that the exclusion of cross-mode con- 
tent makes this representation of generalized matches less 
closely related to a matched-filtered SNR calculation or 
even, potentially, to the maximum likelihood estimator, 
in that it does not account for circumstances where the 
maximum likelihood occurs at an incorrect value for the 
parameters. In such cases, this method will instead ig- 
nore all maxima except the local maximum determined 
from the largest match of like-modes. In most cases, the 
contribution of cross-mode terms to any match calcula- 
tion will be negligible compared to like-mode contribu- 
tions, so this local maximum will also be the global maxi- 
mum. In that case, the optimized parameter choice found 
from this procedure will be consistent with the maximum 
likelihood value, and the resulting match will be a true 
representation of the fraction of recoverable SNR from 
matched filtering. 
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